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APPARATUS AND METHOD FOR SIMULATING THE RECEIVING 
CHARACTERISTIC OF RADIO WAVES 

Background of the Invention 
5 Field of the Invention 

The present invention relates to an apparatus 
for simulating the receiving characteristic of an 
object that receives radio waves in the analysis of 
radio waves transmitted from a radio wave generation 
10 source and a method thereof. 

Description of the Related Art 

Receiving sensitivity based on the directivity 
characteristic of an antenna is the major factor of 

15 a product that receives radio waves from outside, such 
as a cellular phone, a car antenna, etc. These products 
are assumed to be called "equipment under test (EUT)'''. 
Software programs for virtually modeling and 
calculating the directivity characteristic against 

20 radio waves from outside of an EUT are conventionally 
used. 

Fig. lA shows the relationship between a 
transmitting antenna, which is a radio wave generation 
source (wave source) and a receiving antenna included 
25 in an EUT, in such a model. In Fig. lA, the directivity 
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characteristic of a receiving antenna 12 can be checked 
by rotating a transmitting antenna 11 by arbitrary angle 
8 using the receiving antenna 12 as the center and 
calculating voltage Vin at the arbitrary point of the 
5 receiving antenna 12 when electric field is applied 
to the receiving antenna 12 from the transmitting 
antenna 11. This Vin is called the ^^receiving 
sensitivity" of the receiving antenna 12 . In this case, 
the value of Vin varies depending on rotation angle 

10 9, and the directivity characteristic shown in Fig. 
IB is detected. 

In Fig. IB, a coordinate axis 13 indicates the 
direction in which the receiving sensitivity of the 
receiving antenna 12 is maximum, and length from origin 

15 0 up to point P, which is on a curved line 15, on a 
straight line 14 obtained by rotating this coordinate 
axis by arbitrary angle 0 indicates Vm against 0. In 
other words, the curved line 15 indicates the value 
of voltage Vin against each value of angle 0- 

20 As one method for calculating the directivity 

characteristic of the receiving antenna 12 by replacing 
the antenna 12 with an EUT in an arbitrary shape. Moment 
method is known. Moment method is one solution key to 
an integral equation led from Maxwell's 

25 electromagnetic wave equation, as disclosed, for 
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example, in an " Electromagnetic Field Intensity 
Calculation Apparatus" (Japan Patent Laid-open 
Application No . 7-234890 ) , and it is a method for 
dividing an object into many small elements and 
5 calculating current that flows through each element. 
If the current that flows through each element is 
obtained, the voltage at an arbitrary point of the 
object can be calculated. 

To calculate the directivity characteristic of 

10 an EUT by Moment method, a system consisting of a 
transmitting antenna and an EUT must be modeled, the 
system must be divided into small elements and current 
that flows through each element of the EUT must be 
calculated. The simultaneous equations of Moment 

15 method can be given as follows. 

[Zij] [Ij] = [Vi] (1) 

In the above equation, [Zij] is a matrix having 
20 mutual impedance Zij between the i-th and j-th elements 
of the system as an element, [Ij] is column vector having 
current Ij that flows through the j-th element as an 
element and [Vi] is column vector having the voltage 
Vi of the i-th element as an element. Of these, [Vi] 
25 is given as the wave-source voltage of the model and 
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[Zij] is calculated based on the model . [Ij] corresponds 
to the unknown of the simultaneous equations. 

Fig. IC is a flowchart showing a conventional 
current calculation process by such Moment method. A 
5 conventional processing apparatus first reads the data 
of the model consisting of the transmitting antenna 
and EUT (step SI) and calculates mutual impedance Zij 
against the initial value of the rotation angle of the 
transmitting antenna (step S2) . 

10 Then, the apparatus generates the coefficient 

matrix [Zij] of equation (1) , performs the LDU 
factorization of the matrix (step S3) and calculates 
current Ij by forward substitution/backward 
substitution (step S4) . LDU factorization means an 

15 operation to convert a coefficient matrix into the 
product of lower triangular matrix L, diagonal matrix 
D and upper triangular matrix U, and forward/backward 
substitution means a calculation method for 
calculating a solution using these matrices. 

20 Then, the processing apparatus judges whether 

the angle of the transmitting antenna should be changed 
(step S5) . If the angle shouldbe changed, the apparatus 
calculates mutual impedance Zij against a subsequent 
angle (step S2) and performs the processes in and after 

25 step S3. If the current calculation at all angles is 
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completed, in step S5 the apparatus stops the change 
of the angle and terminates the process. 

However, the conventional calculation method 
described above has the following problems. 
5 According to the conventional calculationmethod, 

mutual impedance Zij in the left side of equation (1) 
must be calculated, mutual impedance matrix [Z^j] must 
be reorganized and the simultaneous equations must be 
solved every time the antenna angle is changed. In this 

10 case, it takes a long time to calculate even the value 
of one piece of mutual impedance Zij . Therefore, if 
the number of system elements increases, it takes an 
enormous time to calculate all elements of the mutual 
impedance Z^j composing mutual impedance matrix [Zij] . 

15 Furthermore, if mutual impedance matrix [Zij] is 
reorganized for k angles, the calculation time becomes 
k times as long as that. 

Summary of the Invention 

20 It is an obj ect of the present invention to provide 

a simulation apparatus for improving the simulation 
speed of the receiving characteristic of an object in 
an arbitrary shape that receives radio waves and a 
method thereof. 

25 The simulation apparatus of the present 
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invention comprises first and second current 
calculation devices, a current storage device and an 
output device. The apparatus simulates the receiving 
characteristic of an object that receives radio waves 
5 transmitted from a wave source. 

The first current calculation device calculates 
the current values of the wave source using the 
simultaneous equations of the wave source, which have 
currents that flow through respective elements as 

10 unknowns when the wave source is divided into a 
plurality of elements. The current storage device 
stores the current values of the wave source . The second 
current calculation device calculates the current 
values of the object using the simultaneous equations 

15 of the object, which have currents that flow through 
respective elements and the current values stored in 
the current storage device as unknowns and constants, 
respectively, when the object is divided into a 
plurality of elements and the positional relationship 

20 between the wave source and object changes. The output 
device calculates the receiving characteristic of the 
object based on the current values of the object, and 
outputs the receiving characteristic. 
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Brief Descriptions of the Drawings 

Fig. lA shows transmitting and receiving 
antennas . 

Fig. IB shows the directivity characteristic of 
5 the receiving antenna. 

Fig. IC is a flowchart showing a conventional 
current calculation process. 

Fig. 2A shows the basic configuration of the 
simulation apparatus of the present invention. 
10 Fig. 2B shows an analysis model. 

Fig. 3 shows the simultaneous equations of the 

model . 

Fig. 4 shows the simultaneous equations of the 
current of the transmitting antenna. 
15 Fig. 5 shows the simultaneous equations of the 

current of an BUT. 

Fig. 6 shows the configuration of the simulation 
apparatus . 

Fig . 7 is a flowchart showing the first simulation 
20 process. 

Fig. 8 is a flowchart showing the second 
simulation process. 

Fig. 9 shows the configuration of an information 
processing device. 
25 Fig. 10 shows storage media. 



Descriptions of the Preferred Embodiments 

The detailedpref erred embodiments are described 
below with reference to the drawings. 
5 Fig. 2A shows the basic configuration of the 

simulation apparatus of the present invention. The 
simulation apparatus shown in Fig . 2 A comprises current 
calculation devices 21 and 22, a current storage device 
23 and an output device 24. The apparatus simulates 

10 the receiving characteristic of an obj ect that receives 
radio waves transmitted from a wave source. 

The current calculation device 21 calculates the 
current values of the wave source using the simultaneous 
equations of the wave source, which have currents that 

15 flow through respective elements as unknowns when the 
wave source is divided into a plurality of elements. 
The current storage device 23 stores the current values 
of the wave source. The current calculation device 22 
calculates the current values of the object using the 

20 simultaneous equations of the object, which have 
currents that flow through respective elements and the 
current values stored in the current storage device 
23 as unknowns and constants, respectively, when the 
object is divided into a plurality of elements and the 

25 positional relationship between the wave source and 
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object changes. The output device 24 calculates the 
receiving characteristic of the object based on the 
current values of the object, and outputs the receiving 
characteristic . 
5 The wave source, for example, corresponds to a 

transmitting antenna, and the object, for example, 
corresponds to the receiving antenna or EUT. The 
simulation apparatus generates simultaneous equations 
concerning the wave source and simultaneous equations 

10 concerning the object separately, and calculates 
current values. 

First, the current calculation device 21 
calculates the current value of the wave source by 
solving the simultaneous equations concerning the 

15 currents of a plurality of elements composing the wave 
source, and stores the values in the current storage 
device 23. Then, when the relative-position 
relationship between the wave source and obj ect changes, 
the current calculation device 22 extracts the current 

20 values stored in the current storage device 23, and 
generates simultaneous equations concerning the 
currents of the plurality of elements composing the 
object using those values as constants. Then, the unit 
22 calculates the current values of the object by 

25 solving the simultaneous equations and outputs the 
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values to the output device 24. The output device 24 
calculates the receiving characteristic of the object, 
such as receiving sensitivity, etc., using the received 
current values and outputs the characteristic. 
5 According to such a simulation apparatus, the 

simultaneous equations of a wave source and 
simultaneous equations of an object can be separated 
by regarding the current values of the wave source to 
be constant, and the current values of the wave source 

10 and the current values of the object can be separately 
calculated. In this case, the coefficient matrix of 
the simultaneous equations of the wave source can be 
composed of only mutual impedance between the elements 
of the wave source, and the coefficient matrix of the 

15 simultaneous equations of the object can be composed 
of only mutual impedance between the elements of the 
object . 

Since these pieces of mutual impedance do not 
change even if the relative position of the wave source 

20 against the obj ect changes, one time of the calculation 
of the coefficient matrix is sufficient. In this way, 
there is no need to repeat the calculation of the 
coefficient matrix for each angle of the transmitting 
antenna as in the conventional method, and as a result, 

25 process time can be greatly reduced- 
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For example, the current calculation devices 21 
and 22 shown in Fig. 2A correspond to the current 
calculation unit 53 shown in Fig. 6, which is described 
later, the current storage device 23 shown in Fig. 2A 
5 corresponds to the current storage unit 63 shown in 
Fig. 6, and the output device 24 shown in Fig. 2A 
corresponds to the voltage calculation unit 55 shown 
in Fig. 6. 

Fig. 2B shows the analysis model of a system 

10 consisting of a transmitting antenna and an EUT . In 
the model of Fig. 2B, a transmitting antenna 31 
corresponds to the wave source of radio waves, and it 
transmits radio waves by applying electric field to 
the EUT 32. The EUT 32 corresponds to a car equipped 

15 with a glass antenna, and it receives radio waves 
transmitted from the transmitting antenna 31. In this 
case, equation (1) canbe replaced with the simultaneous 
equations shown in Fig. 3. 

In the coefficient matrix shown in Fig. 3, a 

20 submatrix 41 has mutual impedance ZEi,j ( i=l , n, j = 1 , 

n) among n elements composing the EUT 32 as an element, 
and submatrix 44 has mutual impedance ZAi,j (i=l,..., m, 
j = l,..., m) among m elements composing the transmitting 
antenna 31 as an element. 

25 A submatrix 42 has mutual impedance ZTi,j (i=l,..,, 
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n, j = l,..., m) between an element composing the EUT 32 
and an element composing the transmitting antenna 31, 
as an element, and submatrix 43 has mutual impedance 
ZTi,j (i = l,..., m, j = l,..., n) as an element. 
5 Current lEj (j=l, n) indicates a current that 

flows through each element of the EUT 32, and current 
lAj (j = l, m) indicates a current that flows through 
each element of the transmitting antenna 31. Voltage 
V is the transmitting voltage of the transmitting 

10 antenna 31. 

If there is a sufficient distance between the 
transmitting antenna 31 and EUT32, it is considered 
that the influence on the transmitting antenna 31 of 
current that flows through the EUT 32 is very small. 

15 Therefore, even i f the angle of the transmitting antenna 
31 against the EUT 32 changes, current that flows 
through the transmitting antenna 31 hardly changes and 
can be regarded to be constant. In this case, the 
transmitting antenna 31 can be used as a constant 

20 current source in the calculation of the current of 
the EUT 32. 

Therefore, first, only the transmitting antenna 
31 is modeled, and the currents lAi through lAm of the 
transmitting antenna 31 are calculated. Simultaneous 
25 equations having only the currents lAi through lAm as 
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unknowns can be generated as shown in Fig. 4 using the 
submatrix 44 shown in Fig. 3. Current values obtained 
by solving these simultaneous equations are input as 
known values to simultaneous equations having the small 
5 matrices 41 and 43 as coefficient matrices, and 
simultaneous equations in which only currents lEi 
through lEn of the EUT 32 are unknown, are generated. 
Simultaneous equations concerning currents lEi through 
lEn are generated as shown in Fig. 5. 

10 In Fig. 5, voltage terms on the right side are 

given by the product of current lAj and mutual impedance 
ZTi,j and they vary depending on the angle of the 
transmitting antenna 31. However, since mutual 
impedance matrix [ZEi,j] on the left side does not vary 

15 depending on the angle, one time of this matrix 
calculation is sufficient. 

For example, if it is assumed that n is 
approximately 30, 000-40, 000 and mis approximately 10 0, 
the calculation of mutual impedance ZEi,j requires far 

20 longer time than the calculation of mutual impedance 
ZTi,j . Therefore, by omitting the calculation of mutual 
impedance ZEi,j when the angle is changed, the speed 
of current calculation can be greatly improved. Such 
a simplified calculation method can be generally used 

25 when there is a sufficient distance between a 



transmitting unit, which is a wave source, and a 
receiving unit that receives radio waves. 

Fig. 6 shows the configuration of the simulation 
apparatus based on such a current calculation method. 

The simulation apparatus shown in Fig. 6 
comprises an impedance calculation unit 51, an LDU 
factorization unit 52, a current calculation unit 53, 
a voltage term calculation unit 54, a voltage 
calculation unit 55, a matrix storage unit 61, an 
impedance storage unit 62, a current storage unit 63 
and a voltage term storage unit 64. 

The impedance calculation unit 51 calculates the 
mutual impedance of a given model and stores the 
impedance in the impedance storage unit 62. The LDU 
factorization unit 52 generates a mutual impedance 
matrix having the calculated mutual impedance as an 
element, performs the LDU factorization of the matrix 
and stores the factorized matrix data in the matrix 
storage unit 61. 

The current calculation unit 53 calculates 
currents using necessary data out of data from the 
impedance calculation unit 51, data from the LDU 
factorization unit 52, data from the matrix storage 
unit 61 and data from the voltage term storage unit 
64, and stores the currents in the current storage unit 
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63 . 

The voltage term calculation unit 54 calculates 
the voltage terms shown in Fig. 5 using both the data 
from the impedance calculation unit 51 and data from 
5 the current storage unit 63 , and store the voltage terms 
in the voltage term storage unit 64. 

The voltage calculation unit 55 calculates 
voltage in the prescribed position of the EUT using 
the data from the current calculation unit 53, and 
10 outputs the value as a simulation result. 

Fig. 7 is a flowchart showing the simulation 
process of the simulation apparatus shown in Fig. 6. 

The simulation apparatus first judges whether 
a calculation method that regards the current of the 
15 transmitting antenna to be constant is applicable to 
the given model (step Sll) . For example, the apparatus 
checks whether a distance between the transmitting 
antenna and EUT is equal to or longer than a prescribed 
threshold value. If the distance is equal to or longer 
20 than the threshold value, the apparatus judges that 
this calculation method is applicable. If the distance 
is shorter than the threshold value, the unit judges 
that this calculation method is not applicable. 

If this calculation method is applicable, the 
25 simulation apparatus reads the data of the model 
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consisting of the transmitting antenna and EUT (step 
S12) . The impedance calculation unit 51 calculates the 
mutual impedance ZAi,j of the transmitting antenna only 
and outputs the calculation result to the current 
5 calculation unit 53 (step S13) . Then, the current 
calculation unit 53 generates the simultaneous 
equations shown in Fig. 4 using the received mutual 
impedance ZAi,j, calculates the current lAj of the 
transmitting antenna and stores the current in the 

10 current storage unit 63 (step S14) . 

Then, the impedance calculation unit 51 
calculates the mutual impedance ZEi,j of the EUT only, 
and outputs the calculation result to the LDU 
factorization unit 52 (step S15) . Then, the LDU 

15 factorization unit 52 generates mutual impedance 
matrix [ZEi,j] using the receivedmutual impedance ZEi, 3 , 
performs the LDU factorization of the matrix and stores 
the factorization result in the matrix storage unit 
61 (step S16) . 

20 Then, the impedance calculation unit 51 

calculates the mutual impedance ZTi,j between the 
transmitting antenna and EUT, and outputs the 
calculation result to the voltage term calculation unit 
54. Then, the voltage term calculation unit 54 

25 calculates the voltage terms shown in Fig. 5 using both 
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the received mutual impedance ZTi,j and the current lAj 
stored in the current storage unit 63, and stores the 
voltage terms in the voltage term storage unit 64 (step 
S17) . 

5 Then, the simulation apparatus judges whether 

the angle of the transmitting antenna should be changed 
(step S18). If the angle should be changed, the 
impedance calculation unit 51 calculates new mutual 
impedance ZTi,j against a subsequent angle. The voltage 

10 term calculation unit 54 calculates new voltage terms 
using both the new mutual impedance ZTi,-, and the current 
lAj stored in the current storage unit 63, and stores 
the voltage terms in the voltage term storage unit 64 
(step S17) . Then, when the calculation of voltage terms 

15 at all angles is completed, in step S18, the simulation 
apparatus stops the change of the angle. 

Then, the current calculation unit 53 generates 
the simultaneous equations shown in Fig. 5 using both 
the factorization result of the mutual impedance matrix 

20 [ZTi,3] stored in the matrix storage unit 61 and the 
voltage terms stored in the voltage term storage unit 
64 . Then, the unit 53 calculates EUT current lEj against 
each angle by forward/backward substitution and 
outputs the current to the voltage calculation unit 

25 55 (step S19) . 
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Then, the voltage calculation unit 55 calculates 
EUT voltage against each angle using the received 
current IE-, and outputs the calculation result as the 
receiving characteristic of the EUT (stepS20) . In this 
5 case, the calculation result of the voltage values is 
displayed on the screen, for example, in the form of 
a directivity characteristic graph, as shown in Fig. 
IB. 

If in step Sll the calculation method in which 

10 the current of the transmitting antenna is constant 
is not available, the simulation apparatus calculates 
currents according to the process shown in Fig, IC (step 
S21) and performs the process in step S20. 

According to such a simulation process, one time 

15 of the calculation of EUT mutual impedance ZEi,j, which 
takes the longest time, is sufficient. Therefore, 
process time can be greatly reduced. Since one time 
of the LDU factorization of a mutual impedance matrix 
[ZEi,j] is also sufficient, process speed can be further 

20 improved. 

The effects of this simulation process are 
describedbelow using as an example an analysis in which 
the directivity characteristic of a glass antenna of 
a car is calculated by irradiating radio waves to the 

25 car from a transmitting antenna. If there are 72 
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transmitting antenna angles to be simulated, according 
to the conventional calculationitiethod, first the total 
analysis time is calculated as follows. 

5 Analysis time = (Analysis time per angle) x72 (2) 

According to the calculation method shown in Fig . 
7, the total analysis time is calculated as follow. 

10 Analysis time = (Analysis time per angle) xl 
+ (Calculation time of mutual impedance between antenna 
and EUT)x71 (3) 

In these equations, (Calculation time of mutual 
15 impedance between antenna and EUT) « (Analysis time 
per angle) . Therefore, the analysis time of equation 
(3) can be regarded to be almost equal to analysis time 
per angle . In other words, process speed can be improved 
approximately 72 times as fast as that by adopting the 
20 calculation method shown in Fig. 7. 

Although in the simulation process shown in Fig. 
7, EUT directivity characteristic against the change 
in rotation angle of a transmitting antenna is simulated, 
similarly, an EUT receiving characteristic against the 
25 change in relative position of a transmitting antenna 



20 



against an EUT can also be simulated. 

In this case, the simulation apparatus 
calculates mutual impedance ZTi,j corresponding to a 
new position by changing the position of the 
5 transmitting antenna. Then, the apparatus generates 
EUT simultaneous equations against the new position 
using the mutual impedance ZTi,j, current lAj stored 
in the current storage unit 63 and matrix data stored 
in the matrix storage unit 61, and calculates new 

10 current values. 

As described above, the simulation process is 
effective if there is a sufficient distance between 
a transmitting antenna and an EUT. However, process 
speed can also be improved without such a condition. 

15 For example, the mutual impedance calculation time 
shown in Fig. 3 can be reduced by storing EUT mutual 
impedance ZEi,-, and using the impedance in the current 
calculation at each angle. 

Fig. 8 is a flowchart showing such a simulation 

20 process. 

The simulation apparatus first reads the data 
of a model consisting of a transmitting antenna and 
an EUT (step S31) . The impedance calculation unit 51 
calculates the mutual impedance ZAi,-, of the 

25 transmitting antenna and stores the impedance in the 
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impedance storage unit 62 (step S32) . 

Then^ the impedance calculation unit 51 
calculates the mutual impedance ZEi, j of the EUT, stores 
the impedance in the impedance storage unit 62 (step 
5 S33) , calculates the mutual impedance ZTi,j between the 
transmitting antenna and EUT and outputs the 
calculation result to the LDU factorization unit 52. 

In this case, the mutual impedance ZAi,j and ZEi,^ 
are stored in the impedance storage unit 62 as data 
10 independent from the angle of the transmitting antenna, 
and ZTi,-, is outputted to the LDU factorization unit 
52 as data dependant on the angle of the transmitting 
antenna . 

Then, the LDU factorization unit 52 extracts the 
15 mutual impedance ZAj,-, and ZEi,j stored in the impedance 
storage unit 62, and generates the mutual impedance 
matrix shown in Fig. 3 using those pieces of data and 
the mutual impedance ZTi,j received from the impedance 
calculation unit 51 (step S35) . Then, the unit 52 
20 performs the LDU factorization of thematrix and outputs 
the factorization result to the current calculation 
unit 53 (step S36) . 

Then, the current calculation unit 53 generates 
the simultaneous equations shown in Fig. 3 using the 
25 received analysis result of the mutual impedance matrix . 
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Then, the unit 53 calculates both the current lEj of 
the EUT and the current IA-, of the transmitting antenna 
by forward/backward substitution, and outputs the 
current lEj to the voltage calculation unit 55 (step 
5 S37) . 

Then, the simulation apparatus judges whether 
the angle of the transmitting antenna should be changed 
(step S38) . If the angle should be changed, the 
impedance calculation unit 51 calculates new mutual 

10 impedance ZTi,j against a subsequent angle and outputs 
the calculation result to the LDU factorization unit 
52 (step S34) . Then, the LDU factorization unit 52 
generates a mutual impedance matrix using the received 
mutual impedance ZTi,j, and ZAi,j and ZEa,j stored in the 

15 impedance storage unit 62. Then, processes in steps 
S36 and S37 are repeated based on the newly generated 
mutual impedance matrix. 

If in this way, current calculation at all angles 
are completed, in step S38 the simulation apparatus 

20 stops the change of the angle. Then, the voltage 
calculation unit 55 calculates EUT voltage against each 
angle using the received current XEj and outputs the 
voltages as a receiving characteristic (step S39) . 

According to such a simulation process, as in 

25 the process shown in Fig. 7, one time of the calculation 
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of EUT mutual impedance ZEi,;, is sufficient . Therefore, 
process time can be greatly reduced. 

Although in the example of Fig. 2B, a car is used 
as EUT, in this preferred embodiment, an analysis model 
5 canalsobe generated using an arbitrary obj ect instead 
of the car. Although in this preferred embodiment, LDU 
factorization is used as the solution key to 
simultaneous equations, an arbitrary matrix 
factorization method can also be used instead of the 

10 LDU factorization- For example, LU factorization for 
converting a coefficient matrix into the product of 
lower triangle matrix L and upper triangle matrix U 
can also be used. 

The simulation apparatus shown in Fig. 6 can be 

15 configured using, for example, the information 
processing device (computer) shown in Fig. 9. The 
information processing device shown in Fig . 9 comprises 
a CPU (central processing unit} 71, a memory 72, an 
input device 73, an output device 74, an external 

20 storage device 75, a medium drive device 7 6 and a network 
connection device 77, and those are connected to one 
another by a bus 78. 

The memory 72 includes, for example, a ROM 
(read-only memory) , a RAM (random-access memory) , etc., 

25 and it stores both a program and data to be used for 



24 



the process. The CPU 71 performs necessary processes 
by using the memory 72 and executing the program. 

The impedance calculation unit 51, LDU 
factorization unit 52, current calculation unit 53, 
5 voltage term calculation unit 54 and voltage 
calculation unit 55 that are shown in Fig. 6 correspond 
to a software component described by a program and each 
unit is stored in the specific program code segment 
of the memory 72. 

10 The input device 73 is, for example, a keyboard, 

a pointing device, a touch panel, etc., and is used 
for a user to input instructions and information. The 
output device 74 is, for example, a display, a printer, 
a speaker, etc., and is used to output inquiries and 

15 process results to a user. 

The external storage device 75 is, for example, 
a magnetic disk device, an optical disk device, a 
magneto-optical disk device, a tape device, etc. The 
information processing device stores the program and 

20 data described above in this external storage device 
75, and uses them by loading them into the memory 72, 
as requested. The external storage device 75 can also 
be used as the matrix storage unit 61 , impedance storage 
unit 62, current storage unit 63 and voltage term 

25 storage unit 64 . 
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The medium drive device 7 6 drives a portable 
storage medium 79 and accesses the recorded contents. 
For the portable storage medium 79, an arbitrary 
computer-readable storage medium, such as a memory card, 
5 a floppy disk, a CD-ROM ( compact disk read-only memory ) , 
an optical disk, a magneto-optical disk, etc. are used. 
A user stores the program and data in this portable 
storage medium 79, and uses them by loading them into 
the memory 72, as requested. 
10 The network connection device 77 is connected 

to an arbitrary network, such as a LAN (local area 
network), etc., and transmits/receives data 
accompanying communications. The information 
processing device receives the program and data from 
15 another device, such as a server, etc., via the network 
connection device 77, and uses them loading them into 
the memory 72, as requested. 

Fig. 10 shows computer-readable storage media 
for providing the information processing device shown 
20 in Fig. 9 with both a program and data. The program 
and data that are stored in the portable storage medium 
79 or the database 81 of a server 80 are loaded into 
the memory 72. In this case, the server 80 generates 
a propagation signal for propagating the program and 
25 data, and transmits the propagation signal to the 
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information processing device via an arbitrary 
transmission medium on a network. Then, the CPU 71 
performs necessary processes by using the data and 
executing the program. 
5 According to the present invention, in the 

simulationof the receiving characteristic of an object 
in the case where radio waves are transmitted from a 
wave source to an object in an arbitrary shape, the 
redundant calculation of mutual impedance can be 
10 omitted and as a result, the process speed of simulation 
can be improved. 



